Berry phase mediated topological thermoelectric transport in gapped single and bilayer graphene 
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We consider the anomalous thermoelectric transport in gapped single and bilayer graphene where the gap may 
be due to broken inversion symmetry. In the presence of the gap, non-trivial Berry phase effects can be shown to 
mediate a transverse thermoelectric voltage in response to an applied temperature gradient even in the absence of 
a perpendicular magnetic field. This spontaneous, anomalous Nernst effect is nonzero for non-uniform chemical 
potential in the two inequivalent valleys in the graphene band structure. Conversely, the Nernst response can be 
used to create a valley-index polarization between the two transverse sample edges as in the analogous valley 
Hall effect. 
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I. INTRODUCTION 

Interesting electric transport properties of the Dirac quasi- 
particles of graphene have recently come under focus i'' 2 i 3 i 4 i 5 
Graphene, which consists of a single layer of carbon atoms 
in a bi-partite honeycomb lattice, has a unique band struc- 
ture that includes two inequivalent valleys in the first Bril- 
louin zone of the momentum space. In gapless graphene, the 
valleys are characterized by the conduction and the valence 
bands touching each other at zero energy. Near the valleys, 
the appropriate quasiparticle dispersion relations are linear in 
momentum, similar to those of massless Dirac quasiparticles. 
Inter- valley-scattering is generally suppressed 6,7,8 , and will be 
ignored here, due to the absence of a large enough scattering 
wave vector mixing the two Dirac points. 

In some graphene samples, e.g., epitaxial graphene on SiC 
substrate, a sizable gap has been reported 9 in the quasiparti- 
cle spectrum near the Dirac points. Since the gap approaches 
zero with increasing sample thickness, it has been attributed to 
the effects due to the substrate, which is believed to break the 
sublattice (hence, the space inversion (SI)) symmetry. 9 Note 
that a quasiparticle gap in graphene can also result from con- 
fined geometries such as in the graphene nanoribbons.i&i! In 
bilayer graphene, an applied interlayer voltage bias can lead 
to a gap due to broken SI symmetry. 12 In such systems, the 
effects of the momentum space Berry curvature, 13 which is 
sharply peaked at the two valleys, give rise to non-trivial, 
topological, electric transport phenomena*^*^ In this paper, 
we discuss the effects of the topological Berry phases on 
another class of transport coefficients, the thermoelectric re- 
sponse coefficients. We will focus, in particular, on the topo- 
logical Nernst effect of gapped graphene. In light of the recent 
successful measurements of the regular thermoelectric coeffi- 
cients in graphene j 16 i 17 ' 18 the anomalous, topological Nernst 
response we discuss here acquires particular relevance and 
importance. Although the specific calculations in this paper 
apply to graphene samples with broken SI symmetry, they 
should, more generally, be qualitatively valid in the presence 
of a gap of any origin near the Dirac points. Note that an 
energy gap is ultimately essential for use of graphene as an 
electronic material. 

We calculate the Berry phase supported, spontaneous {i.e., 



it exists even without an external magnetic field) Nernst effect 
of the low energy quasiparticles in graphene, and show that 
it is non-zero for a non-zero valley polarization in the sam- 
ple. While a non-zero valley polarization can be artificially 
created^ it can also be induced by an applied magnetic field 
itself in a generic Nernst measurement set up, 1 - The regular 
Nernst effect in gapped graphene may, therefore, have a sub- 
stantial anomalous contribution because of the induced valley 
polarization. Even in the absence of any valley polarization, 
the spontaneous Nernst signals are non-zero from the individ- 
ual valleys, but are equal and opposite to each other so the 
total signal is zero. Therefore, they can be used to create a 
non-zero valley polarization on the transverse edges, which 
can be technologically important. 

In recent experiments, thermoelectric coefficients, in partic- 
ular, the Nernst coefficient, have been measured for non-zero 
transverse magnetic fields in gapless graphene , 16 > 17 i 18 The 
topological Nernst effect we discuss here exists even in the 
absence of any external magnetic field; it is solely driven by 
the effective magnetic field, given by the Berry curvature, near 
the valley centers in the momentum space. In the thermoelec- 
tric experiments, a temperature gradient — VT, applied along, 
say, the x direction produces a measurable transverse electric 
field. The total charge current in the presence of the electric 
field E and — VT is given by Ji = (TijEj+aij(—djT), where 
<7jj and aij are the electric and the thermoelectric conductiv- 
ity tensors, respectively. Setting J to zero, the Nernst signal, 
defined as e^r = E y J |VT| = pa xy — S xx t&n9H, is mea- 
sured, where a xy is the Nernst conductivity defined via the 
relation J x = a xy (—d y T) in the absence of the electric field, 
p = 1/<J XX is the longitudinal resistance, S xx — E x / |VT| = 
pa xx is the thermopower, and tan## = a xy /a xx is the Hall 
angle. If the second term is small, which is sometimes the 
case,— pa xy completely defines the Nernst signal, but in the 
most general case one should extract a xy from the experimen- 
tal data to compare with our calculations. 

In this paper, the Berry phase mediated topological contri- 
bution to a xy is calculated through two different approaches. 
In the first, path integral, method, we first compute the 
frequency-dependent off-diagonal component of the electric 
conductivity tensor, (7^(0), for the Dirac quasiparticles in 
graphene. At low temperatures, the Nernst conductivity a xy 
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can be related to the zero temperature DC Hall conductivity 
a xy through the Mott relation, 21 Eq. (l22l >. This method has 
the advantage that it produces, as a by-product, the imaginary 
part of the AC Hall conductivity, <r" y , which can be related 
to the possible polar Kerr effec t 22 ' 23 in graphene. Indeed, in 
the case of a non-zero valley-index polarization in the sample, 
the Berry phase induced polar Kerr effect will be non-zero 
and observable. Another advantage of this method is that it 
can give a neat analytical expression of the Nernst conduc- 
tivity at low temperatures. The effects of quasiparticle non- 
linear dispersion relations on the topological transport proper- 
ties of graphene have also been discussed in this framework. 
This method, however, is approximate, and gives the values 
of the Nernst coefficient only at low enough temperatures. In 
a more general method, which is valid also at higher temper- 
atures, the anomalous contribution to a xy can be calculated 
from the coefficient a xy , which determines the transverse heat 
current J h in response to an electric field E: j' x l = a xy E y . 
The coefficients are related by one of the Onsager relations: 
a xy = Ta xy M& This way of computing a xy makes use of 
the Berry curvature in the graphene band structure (see Eq. 
(l24l i) more directly. It is a much more general method, which 
we have used in our plots for the Nernst coefficient as a func- 
tion of temperature in Figs. 2 and 3. Using this method, we 
have been able to go beyond the linear temperature depen- 
dence arising out of the Mott relation. 

The paper is organized as follows: Sec. II gives the Dirac- 
like Hamiltonian for single layer graphene. Section III is 
devoted to the topological Hall conductivity of single layer 
graphene. We calculate the AC Hall conductivity (7^(0) us- 
ing the path integral approach and discuss the possible polar 
Kerr effect for single layer graphene. In the limit £1 = 0, 
we obtain the DC Hall conductivity. The effects of nonlinear 
quasiparticle dispersion relation on the Hall conductivity are 
also discussed. In Sec. V, we calculate the topological Nernst 
conductivity a xy via both the Mott relation and the transverse 
heat current coefficient. We do this for both single and bilayer 
graphene. Finally, Section VI is devoted to the conclusion. 



with 5i,2 = § (l,±V3), S 3 = a(-l,0)x, a is the lattice 
constant. The Hamiltonian (HJ operates on the two-component 
spinor i!f p = (ca p , Cb p ), where ca p , cb p denote the annihila- 
tion operators on the two sublattices A, B. The most interest- 
ing physics occurs around two valleys located at the Brillouin 
zone corners Ki, 2 = 5^ (0, r»), where t z = ±1 labels the 
two valleys. Denoting k = p — K 12 to measure the momen- 
tum from the valley centers, we can rewrite the Hamiltonian 
(D as 



H = 4-0* 
2 



■ Ci (k) a x + C2 (k) a y . 



Here 



Ci(k) = - 
C 2 (k) = t 



2 cos 



2 sin ^ cos 



+ cos (k x a) 



sin (k x a) 



(3) 

(4) 
(5) 



The spectrum of the Hamiltonian consists of two branches 
with the eigenenergies E±(k) = ±A(k), where 



A(k) = 

e = 



2tt 



l + 4cos 2 9 
k v a\/3 



4 cos cos 



3 A/ 2; cl 



A 2 V 
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Close to the valley centers, we can expand d (k) and (2 (k) 
to the first order in k, and rewrite the Hamiltonian as 



H = —o z + — \J z k y a x 



k X <Jy) 



A (k) • a, (7) 



with Ai = 3atr z k y /2, A 2 = -3atk x /2, A 3 = A/2. The 
energy spectrum becomes 



A(k) = |A(k) 



VA 2 + 9aH 2 k 2 



(8) 



II. DIRAC-LIKE HAMILTONIAN 

In the tight-binding approximation, graphene single layer 
with staggered sublattice potential (which breaks the space in- 
version symmetry) can be modeled with the nearest-neighbor 
hopping energy t and a site energy difference A between the 
sublattices 



H 



A/2 C(P) 
C(p) -A/2/' 



(1) 



where 

C(p) 



-te 



ip x a/2 



2cos ^ Ue-*"/ 2 
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III. HALL CONDUCTIVITY 

A. AC Hall conductivity 

Introducing the Lagrangian density, C = iuj n I — H(k), 
where uj n is a fermionic Matsubara frequency, we obtain the 
electron Green's function (for a single valley) as a function of 

k = [iio n , k x , ky ) 



G (k) = £- 



iuj n I + A(k) • ct 



(9) 



Here g{k) — (itUn) 2 — |A(k)| 2 . In order to calculate the 
AC Hall conductivity, we consider the electromagnetic po- 
tential A = (Aq, A x , A y ), where Ao is the scalar potential, 
and A = (A x , A y ) is the vector potential. Expanding C to 
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the first order in A, and integrating out the fermions, we iden- 
tify the coefficient of the Chern-Simons-like term in the effec- 
tive action as the anomalous AC Hall conductivity for a single 
valley^ 



T(k) 



g(un, k) g(uj n + 



T(k) = -2A ■ 



dA dA 



(10) 



(11) 



Here the electromagnetic potential A (q) was assumed to be a 
function of the Fourier variables q — (iSl m , q x ,q y ), with f2 m 
a bosonic Matsubara frequency, and we took the limit q — > 
to arrive at Eq. ( TTOb . Substituting the expression for A(k) 
from the Hamiltonian (O into the formula (fTTI) . we find 



T(k) = 



-9r 2 a 2 i 2 A 



(12) 



Summing over the fermionic Matsubara frequency w„ in 
Eq. ( [Tol l and analytically continuing to the real bosonic fre- 
quency, iSl m — > £1 + i(5, where S is a positive infinitesimal, we 
find, 



a xy (Sl) 



d 2 k T(k) {/ [£ + (k) - M ] - / [£_(k) - M ]} 
(2tt) 2 A(k) [fi + itf - 2A(k)] [fi + iS + 2A(k)] ' 

(13) 

Here — a«) is the Fermi occupation function at temper- 
ature T and /i is the chemical potential. The minus sign be- 
tween the two Fermi functions is due to the opposite sign of 
the Berry curvature for the conduction and the valence bands 
of electrons. Since A(k) = VA 2 + 9a 2 t 2 k 2 /2 and T(k) are 
constants, we can define e = fc 2 /2 and change the integration 
to 



a xy (Sl) 



de 



TFQj) 



2tt E (e) [Sl + iS- 2E (e)] [SI + iS + 2E (e)] ' 

(14) 

where E (e) = ±VA 2 + 18a 2 t 2 e, F (fi) = {/ [E (e) -fi]- 
f[-E 

It is clear that if the two valleys have the same chemi- 
cal potential, [Ik 1 = ^k 2 , the total AC Hall conductivity 
a xy (Sl) = cr^ \st) + <J x ^{Sl) = 0, because of the opposite 
signs of the function T for the two valleys. However, the Hall 
current for an individual valley is nonzero. 

The imaginary part of the AC Hall conductivity, a" y (fl), 
comes from the pole in Eq. ( fT3l at SI = 2E (e). This pole 
represents vertical interband transitions induced by the elec- 
tromagnetic wave. Taking the imaginary part of (1131 1. we find 



TiSPh 



deS[Sl - 2E (e)]F (fj.) . (15) 



Here we have already taken into account the spin degeneracy. 
At T = 0, F (fj) = if E (e) < fi and FQu) = 1 for 
E (e) > /i. Therefore, if SI < 2[i or SI < A, a'^Sl) = 0. If 
SI > 2/i and SI > A, we have 



e 2 A 
2htt' 



(16) 



At finite temperatures, we have a xy (Sl) = jf^F (jjl) if SI > 
A. For unequal population of the quasiparticles in the two 
valleys, the imaginary part of the AC Hall conductivity is non- 
zero. Since the polar Kerr angle, 9k, which is the angle of ro- 
tation of the plane of polarization of a linearly polarized light 
on normal reflection from the sample, is directly proportional 
to cr" y (r2) 22 i 23 , in this case, there will be a non-zero polar Kerr 
effect in graphene. 



B. DC Hall conductivity 

Taking SI = in Eq. ( fT3l . we have for the DC Hall con- 
ductivity, 



Notice that 

n 



cr xy (0) 



T 



d 2 k TF(fx) 
(^O 2 4A 3 (k)' 



9a 2 t 2 Ar z 



4A 3 (k) 2 (A 2 + 9a 2 t 2 k 2 ) 



3/2 



(17) 



(18) 



is just the Berry curvature for an individual valley 15 . We see 
that the Berry curvatures for the two valleys differ only by a 
sign. At zero temperature and At = 0, <J xy (0) — —t z /Att. 
We should multiply the result by two to take account of spin 
degeneracy. Restoring e and %, we have (for a single valley) 

a xy {0) = -e 2 r z /h. (19) 

For nonzero At > A/2, i.e., /i lying in the conductance band, 

<r xy (Q) = -r z e 2 A/2hfi. (20) 

As At becomes very big, 1X^(0) again becomes zero. 

C. Nonlinear spectrum 

The form of the Berry curvature ( fT8b n = — j^j also ap- 
plies to the general Hamiltonian (O without the linear approx- 
imation (i.e., the Hamiltonian (0). In this case, we should use 
A determined by Eq. ©. It is easy to show that 



T (k) = - 



V3a 2 t 2 Ar z 



2c os(| + ^r z )cos(%^) 
+ cos(| - \/3k y ar z ) 



(21) 

The lowest order in k (i.e., k = 0) of OH yields Eq. (fT2l . In 
the presence of the nonlinear part of the spectrum, the sym- 
metry relation between the Berry curvatures in the two val- 
leys is given by (k) = —Hk 2 (~ k), and the total Nernst 
conductivity will be non-zero only for non-zero valley polar- 
ization. In Fig. QJ, we plot the Berry curvature around the 
valley center K2 (t z = —1). We see that the Berry curvature 
is very small for a large k. Therefore, the large k regimes give 
very little contribution to the transport quantity of graphene. 
In fact, for (k x a, k y a) = ggj, the Berry curvature is only 
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FIG. 1: (Color online) (a) Berry curvature II (we plot —II for con- 
venience) around the valley center K%. t = 2.8eV, A = 0.26eV. 
The plot is in the logarithm scale. The dashed line indicates the inte- 
gration region for (b). (b) Plot of the DC conductivity with different 
integration region rCmax k x > k y ^ femax. 

The unit of a xy (0) is 

e 2 /h. 



1% of its peak value. In Fig. QJ), we plot the DC Hall con- 
ductivity at T = and /i = 0, taking account of the nonlinear 
spectrum. Because of the nonlinear spectrum, we cannot take 
the integration region up to A; = oo, the limit that is used 
for the linear spectrum. However, the DC Hall conductiv- 
ity cr xy (Q) at each valley is still very close to —e 2 r z /h (the 
value in Eq. (T% obtained using the linear spectrum) for a 
reasonably large integration region. Therefore the nonlinear 
spectrum does not play an important role on the anomalous 
topological transport properties of graphene, and will be ne- 
glected in the following. 



IV. ANOMALOUS NERNST EFFECT 

A. Single layer graphene 

Since the transverse Hall conductivity can be non-zero for 
gapped graphene, it is natural to ask whether there exists an 
anomalous Nernst effect. For this purpose, at low tempera- 
tures, we can use the Mott relational, 



n 2 kg da xy 
3e d\i 



(22) 



For —A/2 < fj, < A/2, i.e., /i falling in the band gap, a xy is 
independent of [i, and we have -5-^ = 0. Thus, a xy = 0. If 
[i > A/2, we have 



d<r xy _ T z e 2 l 
d[i 2hjj? 



A 



(23) 



We see that 



Similarly, for /x < -A/2, ^ = -'., TW . 
a xy is discontinuous at /i = ±A/2. 

More generally, the anomalous contribution to a xy can be 
calculated from the transverse heat current coefficient a xy . 
This way of computing a xy makes direct use of the Berry 
curvature in the graphene band structure. In the presence of 
the Berry curvature and the electric field, the electron velocity 
acquires the additional anomalous term ?iv k = eE x n(k) 27 . 
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FIG. 2: (Color online) Anomalous Nernst signal for a single valley 
for single layer graphene. The units for /i, T, and a xy are eV, K, and 
kse/Ti, respectively. The parameters t = 2.82eV, A = 0.26eV. (a) 
Contour plot for a xy . (b) Plot of a xy with respect to /i for T = 11.6 
K. Note the discontinuity at /i = A/2, as predicted by the Mott 
relation. The decay for /i > A/2 follows ~ l//i 2 - (c) Plot of a xy 
with respect to /i for T = 360 K. Note the disappearance of the 
discontinuity, (d) Plot of a xy with respect to T for fi — 0.2 eV. 



This velocity multiplied by the entropy density produces the 
transverse heat current coefficient—: 

2e . f d 2 h 
a xy = Ta xy = — V / ——J n„(k) s„(k). (24) 

Here s(k) = -/ k In / k - (1 - / k ) ln(l - / k ) is the entropy 
density of the electron gas, / k = f[E n (k)] is the Fermi dis- 
tribution function, and the sum is taken over both bands. This 
formula for a xy , which we use in the figures|2]and[3] coincides 
with the expression derived from Eq. d22"b only at low temper- 
atures. If = MKai we see again that a xy = 0, since 
the Berry curvature = — Hk 2 - However, the anoma- 
lous Nernst coefficient for each valley remains non-zero. The 
anomalous Nernst effect, for a uniform chemical potential, 
can therefore be used to create equal and opposite valley- 
index polarization on the two transverse edges. This situation 
is similar to the analogous spin Hall effect 29 or possible val- 
ley Hall effect in graphene 15 . For unequal chemical potential 
in the valleys, which can be either artificially created 19 or in- 
duced by a perpendicular magnetic field in a generic Nernst 
set up, there will be a non-zero anomalous contribution to the 
Nernst effect in graphene. The unequal valley chemical poten- 
tial in the presence of an external magnetic field is due to the 
Zeeman-type coupling of the magnetic field with the substan- 
tial magnetic moment associated with the valleys, which can 
be 30 times bigger than the usual spin magnetic moment^. 

Since Eq. d24T i is more general than the Mott relation, and 
may be used to go beyond the linear temperature dependence 
of a xy 



we numerically evaluate Eq. d24b to calculate the 
Nernst coefficient. We plot the anomalous Nernst coefficient 
as a function of /i and T for a single valley in single layer 
graphene in Fig. (f2]). In Fig. [2^, there is a nonzero Nernst 
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FIG. 3: (Color online) Anomalous Nernst signal for a single valley 
for bilayer graphene. The units for /i, T, and a xy are eV, K, and 
kse/h, respectively. The parameters t = 2.82eV, V = 0.2eV, 
tx — OAeV. (a) Contour plot for a xy . (b) Plot of a xy with respect 
to n for T = 11.6 K. Note the dramatic variation of a xy around 
fj, = V/2. (c) Plot of a xy with respect to fi for T = 360 K. (d) Plot 
of a xy with respect to T for /i = 0.2 eV (Solid line), 0. 1 eV (Dashed 
line), 0.091 eV (Dotted line). 



^ = (1 + Tin + 72„ + lln) 

7m = tk/(E n -V/2), 72 „ 

- [t 2 k 2 / {E n - V/2) - {E n - V/2)] /t x , 73 „ 
j 2 tk/ (E n + V/2), and, finally, 



E n 



±- 



2t\ + V 2 + 4t 2 k 2 
±2^A{V 2 + tl)t 2 k 2 + t] 



(27) 



are the eigenenergies. 

In Fig. [3] we plot a xy as a function of [i and T, similar as 
in Fig. |2] Note that, in this case, a xy can change sign in some 
regimes of the parameter space. Such sign changes are due to 
the complicity of the four bands of bilayer graphene. At low 
temperatures, a xy — for fj, < V/2 and changes dramatically 
around fi = V/2. The dramatic change of a xy disappears at 
medium temperatures, but the sign change still remains (see 
Fig- 00- At different chemical potentials, the temperature 
dependence of a xy shows quite different behavior (see Fig. 
[3}l). For uniform chemical potential at the valleys, the total 
anomalous Nernst signal is zero, even though, like in the case 
of single layer graphene, they can be substantial locally in the 
momentum space. 



coefficient even when the chemical potential is inside the band 
gap, while the Mott relation yields a zero Nernst signal in this 
case. In Fig. [2}}, Nernst conductivity a xy increases sharply 
to a maximum, while the Mott relation yields a step jump at 
fi = y . In Fig. [3J;, the step jump is further smoothed out at 
medium temperatures (note that T here is still much smaller 
than A). The deviation of a xy from the linear temperature 
dependence given by the Mott relation is clearly visible in Fig. 



B. Anomalous Nernst effect in bilayer graphene 

The Hamiltonian for a single valley in biased bilayer 
graphene can be written as2& 



H = 
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where V is the interlayer bias voltage, t j_ and t describe inter- 
layer and intralayer nearest neighbor hoppings and </>k is the 
angle of k on the 2D momentum space. The eigenenergies 
and eigenstates can be obtained analytically, which yield the 
Berry curvature for the n-th band 



II, 



ld_ 

k dk 



(26) 



C. Conclusion 

In summary, we discuss the effects of the topological 
Berry phases on the thermoelectric response coefficients of 
the Dirac-like quasiparticles in gapped single and bilayer 
graphene. The gap in the single layer graphene can be due 
to SI symmetry breaking, as claimed, for instance, in epitax- 
ial graphene in Ref. 9 . It may also be due to confined ge- 
ometries, such as in graphene nanoribbon a 10 ' 11 . For bilayer 
graphene, the gap may be due to applied voltage bias 1 —. We 
show that in such samples the anomalous contribution to the 
Nernst effect is non-zero for a non-zero valley polarization. 
We also show that the nonlinear quasiparticle dispersion spec- 
trum does not significantly change the results. A non-zero 
valley polarization may be artificially created^ 9 -, or it can be 
induced by a perpendicular magnetic field due to the substan- 
tial valley magnetic moment 15 in a generic Nernst set up. The 
regular Nernst effect in gapped graphene may, therefore, have 
a substantial anomalous contribution because of the valley po- 
larization. Even in the absence of any valley polarization, the 
spontaneous Nernst effect from an individual valley is locally 
non-zero, opening up the possibility to create a measurable 
valley polarization on the transverse edges of a sample which 
can be technologically very important. 
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